Impurity entanglement entropy in Kondo systems from conformal field theory 
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The entanglement entropy in Kondo impurity systems is studied analytically using conformal 
field theory. From the impurity contribution to the scaling corrections of the entanglement entropy 
we extract information about the screening cloud profile for general non-Fermi-liquid fixed points. 
By also considering the finite-temperature corrections to scaling of the von Neumann entropy we 
point out a direct connection between the large-distance screening cloud profile and thermodynamic 
observables such as the specific heat. 
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Introduction. Entanglement — nonseparability of 
states — lies at the very heart of quantum theory. It has 
been recognized as the crucial resource needed for per- 
forming quantum computing and teleportation, but has 
also entered as an important concept in a wide range of 
fields spanning from black holes to biological systems ii 
Entanglement entropy, a measure of the entanglement 
between two parts of a quantum many-body system, has 
been used as a theoretical tool to study quantum phase 
transitions as well as characterizing topological phases 
of matter— and developing numerical algorithms)^ A par- 
ticularly important result has been for one-dimensional 
(ID) critical systems, where conformal field theory yields 
a universal prediction for the scaling of the entanglement 
entropy4i^ When such a system has a boundary, it pro- 
vides a framework for the description of universal features 
of quantum impurity systemsi^ Important information 
can be obtained by studying the subleading corrections 
to the scaling of the entanglement entropyji suggesting 
a new perspective on the long-standing problem of the 
evasive "Kondo screening cloud" <^ 

Regarding the question of measuring entanglement en- 
tropy experimentally, it needs to be related to observable 
quantities. Recent attempts have focused on the connec- 
tion between entanglement and fiuctuationsj^*^ In par- 
ticular, for critical one-dimensional systems with bound- 
aries there is not only the same logarithmic scaling, but 
also the corrections to scaling of the entanglement en- 
tropy and those of the particle (or spin) fluctuations show 
interesting similarities^^ 

In this paper, we employ boundary conformal field the- 
ory (BCFT) to give a unified picture of the entangle- 
ment entropy in quantum impurity problems. In partic- 
ular, we unveil another type of connection between zero- 
temperature entanglement and fiuctuations, this time 
finite-temperature thermodynamic fluctuations (in the 
form of the impurity speciflc heat), at the level of scal- 
ing corrections. As we shall argue below, this provides 
an experimental inroad to the study of Kondo screen- 
ing clouds. The same framework allows us to extend 
the Fermi-liquid analysis of S0rensen et al^ of the large- 
distance behavior of the entanglement entropy generated 
by a Kondo impurity, to the whole range of non-Fermi- 
liquid flxed points in the large family of Kondo systems. 



For the Fermi-liquid case, there is an intuitive picture 
relating the impurity contribution to the entanglement 
entropy and the form of the Kondo screening cloudi^ In a 
valence-bond basis, the screening cloud can be thought of 
as the distribution of singlet bonds between the impurity 
and conduction electrons, hence the valence bond entan- 
glement cntropyii provides the connection. Now, there is 
no such intuitive picture for the non-Fermi-liquid case of 
the more complex Kondo systems, but the analogy pro- 
vides a rationale for regarding the impurity entanglement 
entropy as a good description of the screening cloud pro- 
flle even when the (partial) screening is carried out by the 
highly non-trivial objects that make up the non-Fermi- 
liquid ground state. Importantly, our BCFT result which 
establishes an exact relation between the impurity con- 
tributions to the entanglement and the specific heat in 
terms of scaling corrections of the von Neumann entropy 
then shows that there is a direct connection between the 
large-distance behavior of the screening cloud and an ex- 
perimental observable. 

von Neumann entropy from BCFT. Entanglement en- 
tropy in a zero-temperature bipartite system is en- 
coded in the von Neumann entropy. For a critical 
ID semi-infinite bipartite system at inverse tempera- 
ture /3, the scaling of the von Neumann entropy S{r) = 
—Trp{r) In p{r) of the reduced density matrix p(r) for an 
interval of length r at the boundary of the system is given 
bv^ 
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Here c is the central charge. In g is the boundary en- 
tropyr^ e is a short-distance cutoff, and c' is a non- 
universal constant. This relates the logarithmic scal- 
ing in r of the zero-temperature entanglement entropy 
to the scaling of the extensive thermodynamic entropy 
when r//3 — > oo. We will now see that there is a some- 
what similar relation also for the scaling corrections of 
S(r), denoted by "..." in Eq. ([T]). These corrections are 
governed by the irrelevant operators (in the language of 
the renormalization group) in the BCFT that describes 
the critical properties of the system, and will have both 
bulk and boundary contributions. For quantum impu- 
rity systems, the fixed-point properties can be described 
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by a BCFT where the impurity has been reduced to a 
specific boundary condition and boundary operator con- 
tent.— Therefore, we see from Eq. ^ that the impurity 
contribution Simp{r) to the von Neumann entropy S{r) 
is 



Simpir) = ln.g + 



(2) 



Now is the boundary contribution to the scahng 
corrections of S{r), which is governed by the irrelevant 
boundary operators {</){,} in the theory. 

The von Neumann entropy is calculated as S{r) = 
lim„^i(l - lnTrp"(r). When perturbing the BCFT 
with an irrelevant boundary operator, H — Hbcft + 
X(j)b, the leading correction to S{r) is generically of sec- 
ond order in the scaling field A, and can be written as 
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Here jdz/dwl comes from the mapping z i— > w that 
takes the n-sheeted Riemann surface TZn in the finite- 
temperature geometry (representing Trp") to the finite- 
temperature strip {w = T + iy \ — /3/2 < t < /3/2, y > 0} 
in the complex upper half-plane The cutoff in the 
integral is given by |ti — T2I > €l\dz j dw\^^^-^. Eq. ([3]), 
where Xb is the scaling dimension of (j)b, follows from the 
analysis in Ref. [isl generalized to finite temperature. 

In the limit r / fi 0, we get Idz/dwj oc r+0{{r / fiY). 
This gives, as expected from the zero-temperature result 
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up to terms 0{(r / (3)^). A marginally irrelevant pertur- 
bation generates a leading correction ^ (In^)"'^, while an 
exactly marginal perturbation describes new fixed points 
and hence changes Simp by a constant. The above correc- 
tions are all second-order in A; the only operator which 
gives a non-vanishing first-order correction is the stress- 
energy tensor which gives a correction^ ^ r^^. 

It is also possible to extract the behavior in the limit 
r//3 — > 00 (still at low temperature, i.e. large /3, but tak- 
ing r > ^). We then get that |dz/dw| oc /3 0(^5"^), 
up to terms 0{e^^'^'^/^). Changing integration variable 
to M = tan(7r|Ti — T2|//3) and expanding the integrand in 
Eq. ([3]) in a power series at the divergence at m = 0, gives 
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where /i(Ti) are regular functions. When Xb > 3/2, we 
must use the short-time cutoff, which in the u variable 



becomes u > Tre/ {/3\dz /dw\w=ri) ~ e//3^- Then the lead- 
ing /3 dependence goes as 5 Simp ^ p^-^^b pi^^-Q _ p-i 
when Xb > 3/2, and SSimp ~ ln/3 = I3~^ In when 

Xb = 3/2. When 1 < Xb < 3/2, we see from Eq. ^ 
that the integral ([3]) converges. Hence the leading /3 de- 
pendence comes from the prefactor /J^"^^*" arising from 
jdz/duij cx p. Summarizing, in the limit r//3 — > cx3 we 

get Simp 



\ng + SSimp, with 
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to 0{f3-^) in f3 and ©(e.-^'^''/'^) in r//3. 

The results in Eq. ([6]) bear a close resemblance to 
the well-known expressions for the impurity specific heat 
Gimp at criticalityJ^iii In fact, as Cimp is related to the 
thermodynamic impurity entropy Sfjjip via the relation 
Cimp = dSf^p/dp, they describe the same power law, 
and one has the leading behavior 



9^'* = In 
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if Xb > 3/2, 
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as /3 — !• 00, where Ai, A^ and A3 are constants. Thus the 
von Neumann and thermodynamic impurity entropies 
have the same form on their leading scaling corrections 
in the limit r/ (3 ^ 00 at low temperature. However, the 
amplitudes of these scaling corrections are different, as 
the von Neumann entropy acquires an additional ampli- 
tude factor from the mapping from the Riemann surface, 
not present in the thermodynamic entropy. 

Impurity entanglement entropy in Kondo systems. At 
zero temperature the von Neumann entropy measures the 
entanglement between the two parts of the bipartite sys- 
tem, and the impurity part Simp{r) in Eq. ([2]) is then 
referred to as the impurity entanglement entropy.— We 
can therefore use the zero-temperature results in Eq. (|4]) 
to predict the impurity entanglement entropy in vari- 
ous Kondo impurity models. They all share the com- 
mon feature that the only significant zero-temperature 
length scale is the Kondo length S^k ^ vf/Tk: the char- 
acteristic length scale at which screening is supposed to 
occur (here is the Fermi velocity and Tk the Kondo 
temperature). In particular, the distance dependence of 
Simp can only come through r/^K- Note that when the 
models describe two- (three-) dimensional systems the 
size r of the block at the boundary will correspond to 
the radius of a disc (sphere) centered at the impurity 
(or the midpoint between the impurities when they are 
two) . Thus the impurity entanglement entropy measures 
the impurity-generated entanglement between the part of 
the system within radius r (including the impurity) and 
the rest of the system. It appears as a natural measure 
of the shape of this screening cloud since it captures the 
spatial distribution of the entanglement from the impu- 
rity. Compared to other ways of probing Kondo screening 
with entanglement)^ the impurity entanglement entropy 



3 



has the advantage of allowing analytieal results based on 
BCFT. 

The BCFT approach,— where the model is reduced 
to one spatial dimension with the impurity as a special 
boundary condition, is only valid in the limit r 3> £,k- 
Hence we expect the BCFT prediction for Simp at zero 
temperature to describe the large-distance decay of the 
Kondo screening cloud. By our exact analysis above, the 
leading r dependence of Simp in Eq. (|4]) has the same 
form as the leading /3 dependence of Cimp, so the large- 
distance zero-temperature profile of the screening cloud 
is encoded by the impurity specific heat. The relation to 
the impurity entanglement entropy is less direct for the 
impurity susceptibility, which is governed by the same 
irrelevant boundary operator but through coupling to the 
bulk spin operator We now illustrate our result for a 
large class of Kondo impurity systems. 

The Kondo model. The original Kondo model de- 
scribes a band of conduction electrons interacting with 
a single magnetic s = 1/2 impurity at the origin: 

Hk^ 5^6(fc)4"^,~„ + J5.^^ta(o)^^^(0). (8) 

At the low-temperature Kondo screening fixed point, 
lug = 0, and the leading irrelevant boundary operator is 
the stress-energy tensor T. As found in Ref. adding the 
boundary perturbation 5H = — (^jf/2)T(0) to the fixed 
point Hamiltonian Hbcft gives Simp{r) = 7r^if/(12?') 
for r ^ , in agreement with DMRG results and con- 
sistent with a valence-bond picture. 

Breaking particle-hole symmetry introduces the ex- 
actly marginal charge current operator as a perturbation 
at the fixed point, only giving a constant shift of Sim.p{r). 

The two -impurity Kondo model. Adding a second spin- 
1/2 impurity to the Kondo model gives the two-impurity 
Kondo model (TIKM) 

Htikm = ^ + [ Sc{ri) ■ Si 

+ Sc{r2)-S2]+KSi-S2, (9) 

where Sc(f) = (1/2) Ea,,3 The BCFT 

solution of the model was found in Ref. [l^ In short, the 
model features an unstable fixed point at K = K^ ~ TV, 
where the system undergoes a quantum phase transition. 
At this fixed point In g = In v^, and the leading irrele- 
vant boundary operator allowed by symmetry to appear 
as a perturbation is £-ie, the Virasoro first descendant 
of the e field, with scaling dimension Xh = 3/2. How- 
ever, being a Virasoro first descendant, L_ie will not give 
any contribution neither to the entanglement entropy nor 
any finite-temperature properties to any order in pertur- 
bation theory. This follows from the evaluation of in- 
tegrals of the type dT{L_ie)n„ = dTdr{e)Ti„ 
which vanish due to the periodicity of the boundary of 
the n-sheeted Riemann surface TZn- The leading cor- 
rection to the impurity entanglement entropy therefore 



comes from the stress-energy tensor T, precisely as in 
the single-impurity case. Thus, adding T as a boundary 
perturbation to the fixed-point Hamiltonian Hbcft, 

H = Hbcft + ay/UL-ie{{))-h^KT{Q), (10) 

with a and h dimensionless constants, gives the impurity 
entanglement entropy in the limit r 3> as 

5™p(r) = In V2 + 7r6a'/(60- (H) 

Breaking the SU(2) spin-rotational symmetry or the 
parity symmetry introduces a boundary operator with 
scaling dimensioni^ xi, = 3/2. Particle- hole symmetry 
breaking can be either relevant or exactly marginal. Si- 
multaneously breaking the particle-hole and parity sym- 
metries adds another Xb = 3/2 boundary operator 

The multichannel Kondo model. The Hamiltonian 
is obtained by adding a channel index i = l,2,...,fc 
to the electrons in Eq. ([8]), so that the spin-s im- 
purity interacts with k degenerate bands of conduc- 
tion electrons. When k > 2s, the system is 
governed by an overscreened non-Fermi-liquid fixed 
point, corresponding to a boundary entropy^! Ing = 
\n[sm{TT{2s + l)/{2 + k))/sm{n/{2 + k))]. The leading 
irrelevant boundary operator has scaling dimensioni^ 
xi, ~ l-t-2/(2-|-A:), in agreement with the original 
Bethe Ansatz solution.— Thus, when k = 2 the distance- 
dependent part of Simpij) falls off as ~ iS.K /^) ln(f/$A')- 
For arbitrary k > 2, where 1 < Xh < 3/2, one gets 

S.mpir) = In .9 + A (r/^K )-*/('+'=), (12) 

for r ^ ^if, with A a constant. We see that the more 
channels that are added, the more long-range entangle- 
ment appears to be generated from the screening of the 
impurity by the composite soliton-like objects formed in 
the non-Fermi liquid.— The screening cloud falls off with 
an anomalous power law, a result rather similar to the 
one obtained by Barzykin and Aflleck^^ when defining it 
as the form of the equal-time spin-spin correlator. 

Channel asymmetry is a relevant perturbation,— 
whereas particle-hole symmetry breaking is exactly 
marginal The effect of spin-rotational symmetry break- 
ing depends on s and fcj^ however for the special cases 
where the perturbation is irrelevant the only difference is 
a change in the constant A in Eq. (fT2|) 

Generalizing the multichannel Kondo model by extend- 
ing the spin symmetry group from SU (2) to SU (N) gives 
the multichannel SU(N) Kondo modelM The leading ir- 
relevant operator at the overscreened non-Fermi-liquid 
fixed point now has scaling dimension Xb = l+N/{N+k). 
Hence 5™p(r) = In.g + A{r/^K)-^^^^''+''^ for r » Ck, 
with A a constant and In g depending on the particular 
representation of SU{N). 

The two-impurity, two-channel Kondo model features 
a continuous family of non-Fermi-liquid fixed points,— 
allowing for scaling-correction terms in Simpi^) from 
boundary operators with scaling dimensions ranging all 
the way from a;b = 1 to infinity. 
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The final case we consider is a Kondo impurity in a 
Luttinger liquid. At zero temperature \ng = 0,— but the 
leading perturbation typically has^^ 1 < < 3/2 and 
Eq. (HI) then predicts that S.^^pir) ~ 

Discussion. We have provided a unified picture for the 
impurity entanglement entropy Simp for the large class 
of Kondo models in terms of the scaling corrections aris- 
ing from the boundary operators of the corresponding 
BCFT, vahd for r ^k- When r//3 — !• oo the scaling 
part of the von Neumann entropy approaches the thermo- 
dynamic entropy?^ and it has been argued that this is ex- 
pected on general grounds 4 Our analysis shows that this 
reasoning can also be applied to the scaling corrections, 
but only in the following precise sense: The impurity von 
Neumann entropy and thermodynamic entropy have the 
same leading power-law dependence of /S at low temper- 
atures, however with different amplitudes. This connects 
the exponents for the temperature scaling of the impu- 
rity specific heat with those for the large-distance scaling 



of the zero-temperature impurity entanglement entropy, 
through the von Neumann entropy. 

The result suggests that the large-distance profile of 
a zero-temperature Kondo screening cloud can be read 
off from the impurity specific heat at low temperature, 
an experimentally accessible observable. For the orig- 
inal single- impurity Kondo model, Simp(r) was indeed 
found^ to follow the same power law as the impurity spe- 
cific heat CimpiP), showing that this also holds for the 
first-order correction from the stress-energy tensor. This 
is also what we see when considering the BCFT predic- 
tions for a number of more complex Kondo impurity mod- 
els, thereby demonstrating the close connection between 
Kondo screening, entanglement and thermodynamics in 
quantum impurity systems. 
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